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Abstract 



Measurements from the CMS experiment at the LHC of dihadron correlations for 
charged particles produced in PbPb collisions at a nucleon-nucleon centre-of-mass 
energy of 2.76 TeV are presented. The results are reported as a function of the par- 
ticle transverse momenta (pj) and collision centrality over a broad range in relative 
pseudorapidity (A;/) and the full range of relative azimuthal angle (Acp). The ob- 
served two-dimensional correlation structure in An and Acp is characterised by a nar- 
row peak at (Atj, Acp) (0,0) from jet-like correlations and a long-range structure 
that persists up to at least \An\ = 4. An enhancement of the magnitude of the short- 
range jet peak is observed with increasing centrality especially for particles of pj 
around l-2GeV/c. The long-range azimuthal dihadron correlations are extensively 
studied using a Fourier decomposition analysis. The extracted Fourier coefficients 
are found to factorise into a product of single-particle azimuthal anisotropies up to 
pj ~ 3-3.5 GeV/c for at least one particle from each pair, except for the second-order 
harmonics in the most central PbPb events. Various orders of the single-particle az- 
imuthal anisotropy harmonics are extracted for associated particle pj of 1-3 GeV/ c, as 
a function of the trigger particle pj up to 20 GeV/c and over the full centrality range. 
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1 Introduction 

Measurements of dihadron correlations are a well established technique for studying the prop- 
erties of particle production in the high density medium created in heavy ion collisions. Early 
results from PbPb collisions at the Large Hadron Collider (LHC) (TJ |2| extended these studies 
into a regime of much higher beam energies as compared to those from the Relativistic Heavy 
Ion Collider Il3l-[T0)|. These results complement other LHC measurements of medium proper- 
ties, including a large deficit of charged particles at high-px fTT| and the observations of an 
enhanced fraction of dijets with very asymmetric energies |l2l[T3| . 

The Compact Muon Solenoid (CMS) experiment at the LHC has studied dihadron correla- 
tions over a broad range of relative azimuthal angles (|A</>|) and pseudorapidity (|Aj/|, where 
rj = — ln[tan(0/2)] and 6 is the polar angle relative to the counterclockwise beam axis) in the 
most central PbPb collisions at a nucleon-nucleon centre-of-mass energy (-Js NN ) of 2.76 TeV 
Concentrating on large \At]\, these results established the properties of the so-called "ridge", 
an enhancement of pairs with \A(p\ 0, including in particular a strong dependence of the 
ridge yield on the transverse momentum (pj) of the two hadrons. While a variety of theoretical 
models have been proposed to interpret the ridge phenomena as a consequence of jet-medium 
interactions fl4] - [T9| , recent theoretical developments indicate that, because of event-by-event 
fluctuations in the initial shape of the interacting region, sizable higher-order hydrodynamic 
flow terms could also be induced, e.g., triangular flow |l2"0Tl2"5| . The triangular flow effect will 
contribute to the dihadron correlations in the form of a cos(3A</>) component, which also gives a 
maximum near-side correlation at A<p « 0, similarly to the elliptic flow contribution. It has been 
proposed that by taking into account various higher-order terms, the ridge structure could be 
described entirely by hydrodynamic flow effects |20|. To investigate this possibility, a Fourier 
decomposition of the CMS data at large \Arj\ was performed, again finding a dependence on 
pj 0. Similar results, although with a smaller \At]\ gap, have been reported by ALICE 0. 
The observations by CMS of a ridge-like structure in very high multiplicity proton-proton (pp) 
collisions at a centre-of-mass energy of 7 TeV [26] may also challenge the interpretation of these 
long-range correlations. 

This paper, expanding on previous CMS results [1], presents dihadron correlation measure- 
ments from PbPb collisions at ^/s NN = 2.76 TeV acquired in 2010, for all collision centralities 
and over a broader range of hadron pj. As in Refs. (H|26j, the yield of particles (binned in pj) 
associated with a trigger particle (also binned in pj) is extracted as a function of their relative 
pseudorapidity and azimuthal angle. Such a study of hadron pairs in either the same or dif- 
ferent pj ranges can reveal important information about the production of particles and their 
propagation through the medium. A Fourier decomposition technique is used to quantify the 
long-range azimuthal correlations. The potential connection between the extracted Fourier co- 
efficients from the correlation data and the azimuthal anisotropy harmonics for single particles 
is investigated. This measurement provides a comprehensive examination of the centrality and 
transverse momentum (1 < pj < 20GeV/c) dependencies of the short-range (\At]\ < 1) and 
long-range (2 < \At]\ < 4) dihadron correlations in PbPb collisions at LHC energies, as well as 
the relationship between these two-particle correlations and single-particle angular distribu- 
tions. These results provide extensive input to the interpretation of these observables in terms 
of broad theoretical concepts such as hydrodynamic flow and quantitative models of particle 
production and propagation in the high-density medium. 

The detector, the event selection and the extraction of the correlation functions are described in 
Sections [2] and [3} while the extracted results are described in Sections |4] and [5j 
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3 Data and Analysis 



2 CMS Detector 

The ability of the CMS detector (all components of which are described in Ref. ||27| ) to extract 
the properties of charged particles over a large solid angle is particularly important in the study 
of dihadron correlations. This study is based primarily on data from the inner tracker contained 
within the 3.8 T axial magnetic field of the large superconducting solenoid. The tracker consists 
of silicon pixel and strip detectors. The former includes 1 440 modules arranged in 3 layers, 
while the latter consists of 15148 modules arranged in 10 (11) layers in the barrel (endcap) 
region. The trajectories of charged particles can be reconstructed for pj > lOOMeV/c and 
within \f]\ < 2.5. 

The field volume of the solenoid also contains crystal electromagnetic and brass/ scintillator 
hadron calorimeters. Although not included in the present results, muons are detected using 
gas-ionisation counters embedded in the steel return yoke. In addition to these components in 
and around the barrel and endcap of the solenoid, CMS also has extensive forward calorimetry. 
In the right-handed coordinate system used by CMS, the %-, y~, and z-axes are aligned with 
the radius of the LHC ring, the vertical direction, and the counterclockwise beam direction, 
respectively, with the origin located at the centre of the nominal interaction region. 

For PbPb collisions, the primary minimum-bias trigger uses signals from either the beam scin- 
tillator counters (BSC, 3.23 < \t]\ < 4.65) or the steel/ quartz-fibre Cherenkov forward hadron 
calorimeters (HF, 2.9 < \rj\ < 5.2). Coincident signals from detectors located at both ends of 
the detector (i.e., a pair of BSC or a pair of HF modules) are required. Events due to noise, 
cosmic-ray muons, double-firing triggers, and beam backgrounds are suppressed by further 
requiring the presence of colliding beam bunches. The fraction of inelastic hadronic PbPb col- 
lisions accepted by this primary trigger is (97 ±3)% IIT31 . 

3 Data and Analysis 

The procedure used in the present analysis follows that described in the previous CMS corre- 
lation paper [1J. Offline event selection requires a reconstructed vertex with at least two tracks 
(i.e., at least one pair of charged particles). This vertex must be within 15 cm along the beam 
axis relative to the centre of the nominal collision region and within 0.02 cm in the transverse 
plane relative to the average position of all vertices in a given data sample. In addition, var- 
ious background events (for example beam-gas and beam-halo collisions, cosmic muons, and 
large-impact-parameter electromagnetic collisions) are suppressed by requiring at least three 
signals in the HF calorimeters at both positive and negative rj, with at least 3 GeV of energy in 
each signal. 

The analysis is based on a data sample of PbPb collisions corresponding to an integrated lumi- 
nosity of approximately 3.9 fib -1 p8ll29| , which contains 30 million minimum-bias collisions 
after all event selections are applied. The pp data at yjs = 2.76 TeV, the reference for comparison 
to the PbPb data, were collected during a short low-energy LHC run at the end of March 2011. 
Minimum-bias-triggered pp events corresponding to an integrated luminosity of 520 f/b -1 are 
selected for this analysis. 

The energy released in the collisions is related to the centrality of the heavy ion interactions, i.e., 
the geometrical overlap of the incoming nuclei. The event centrality is defined as the fraction 
of the total cross section, starting at 0% with the most central collisions (i.e., smallest impact 
parameter). This fraction is determined from the distribution of total energy measured in both 
HF calorimeters. 
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The event centrality can be correlated with the total number of nucleons in the two Pb nuclei 
that experienced at least one inelastic collision, N par t- The average values of N par t for the var- 
ious centrality bins used in this analysis are given in Table [T] The N par t values are obtained 
using a Glauber Monte Carlo (MC) simulation I501I3TJ with the same parameters as in Ref. Ifl3| . 
These calculations are translated into reconstructed centrality bins using correlations between 
Npart arid the measured total energy in the HF calorimeters, obtained from fully simulated MC 
events. The systematic uncertainties on the N par t values in Tableware derived from propagation 
of the uncertainties in the parameters of the Glauber model. More details on the determination 
of centrality and N par t can be found in Refs. |T3ll32ll33| . 

Table 1: Average N part values for each PbPb centrality range used in this paper. The values are 
obtained using a Glauber MC simulation with the same parameters as in Ref. IT3"I . 



Centrality 


0-5% 


5-10% 


10-15% 


15-20% 


20-25% 


25-30% 


(Npart) 


381 ±2 


329 ±3 


283 ±3 


240 ±3 


203 ±3 


171 ±3 


Centrality 


30-35% 


35^0% 


40-50% 


50-60% 


60-70% 


70-80% 


(N pa rt) 


142 ±3 


117±3 


86.2 ±2.8 


53.5 ± 2.5 


30.5 ± 1.8 


15.7 ±1.1 



The reconstruction of charged particles in PbPb collisions is based on signals in the silicon pixel 
and strip detectors, similarly to the reconstruction for pp collisions [34|. However, a number 
of settings are adjusted to cope with the challenges presented by the much higher signal den- 
sity in central PbPb collisions. A set of tight quality selections are imposed on the collection 
of fully reconstructed tracks to minimise the contamination from misidentified tracks. These 
include requirements of at least 13 signals on the track, a relative momentum uncertainty of 
less than 5%, a normalized x 2 of less than 0.15 times the number of signals, and transverse 
and longitudinal impact parameters of less than three times the sum in quadrature of the un- 
certainties on the impact parameter and the primary vertex position. Studies with simulated 
MC events show that the combined geometrical acceptance and reconstruction efficiency for 
the primary-track reconstruction reaches about 60% for the 0-5% most central PbPb collisions 
at pj > 2GeV/c over the full CMS tracker acceptance (\t]\ < 2.4) and 66% for \f]\ < 1.0. The 
fraction of misidentified tracks is about 1-2% for \tj\ < 1.0, but increases to 10% at \tj\ « 2.4 
for the 5% most central PbPb collisions. For the peripheral PbPb events (70-80%), the overall 
tracking efficiency improves by up to 5%, with a much lower fraction of misidentified tracks. 

The analysis of dihadron angular correlations in this paper follows exactly the procedure es- 
tablished in Ref. [1J. Any charged particle associated with the primary vertex and in the range 
| rj\ < 2.4 can be used as a trigger particle. A variety of bins of trigger transverse momentum, 
denoted by are considered. There can be more than one such trigger particle in a single 
event and their total multiplicity in a particular data sample is denoted by N tr ig. Within each 
event, every trigger particle is then paired with all of the remaining particles (again within 
\rj\ < 2.4). As for the trigger particles, these associated particles are binned in transverse mo- 
mentum (pj ssoc ). The differential yield of associated particles per trigger particle is given by 

Nag dAtjdAcp v ' B{Ar],A<p) y ' 

where N pair is the total number of correlated hadron pairs. The functions S(At],Acp) and 
B(At],A(p) are called the signal and background distributions, respectively. The value of the 
latter at At] = and A<p = (B(0,0)) is a normalisation factor. 
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3 Data and Analysis 



The signal distribution is the per-trigger-particle yield of pairs found in the same event, 

1 d 2 N same 

where IV same is the number of such pairs within a (At],A(p) bin. The background distribution is 
found using a mixed-event technique, wherein trigger particles from one event are combined 
(mixed) with all of the associated particles from a different event. In the analysis, associated 
particles from 10 randomly chosen events are used. The result is given by 

1 d 2 N mix 

B(Aw,A<*) = -p-=^- — , (3) 

NtngdAT/dA^' v ; 

where N mix denotes the number of mixed-event pairs. This background distribution repre- 
sents the expected correlation if the only effects present were random combinatorics and pair- 
acceptance. 

The value of B ( At], A(p) at At] = and A(p = (with a bin width of 0.3 in At] and 71 / 16 in A<p) is 
used to find the normalisation factor B(0, 0). In this case, the two particles have the maximum 
possible pair acceptance since they are travelling in essentially the same direction. The extent to 
which the background distribution at larger angular separation is smaller than this value (more 
specifically the ratio B(0, 0)/B(At], A(p)) can be used to determine the pair acceptance correc- 
tion factor. Multiplying the signal distribution by this ratio gives the acceptance-corrected per- 
trigger-particle associated yield. Since the distributions should, in principle, be symmetric, the 
statistical precision is maximised by filling only one quadrant using the absolute values of At] 
and A(p. For illustration purposes only (for example, see Fig. [IV, the other three quadrants are 
filled by reflection, giving distributions that are symmetric about (At], Acp) = (0, 0) by construc- 
tion. The pair acceptance decreases rapidly with At] and so, to avoid large fluctuations due to 
statistical limitations, the distributions are truncated at \At]\ =4. The analysis is performed 
in twelve centrality classes of PbPb collisions ranging from the most central 0-5% to the most 
peripheral 70-80%. Within each centrality range, the yield described in Eq. [l]is calculated in 
2 cm wide bins of the vertex position (z vtx ) along the beam direction and then averaged over 
the range |z vtx | < 15 cm. 

When filling the signal and background distributions, each pair is weighted by the product of 
correction factors for the two particles. These factors are the inverse of an efficiency that is a 
function of each particle's pseudorapidity and transverse momentum, 

where A(t],pj) is the geometrical acceptance, E(t],pj) is the reconstruction efficiency, and 
F(t], pj) is the fraction of misidentified tracks. The effect of this weighting factor only changes 
the overall scale but not the shape of the associated yield distribution, which is determined by 
the signal-to-background ratio. 

As described in Ref. [1], the track-weighting procedure is tested using MC events generated 
with HYDJET II35II (version 1.6) propagated through a full detector simulation. The tracking 
efficiencies themselves are checked using simulated tracks embedded into actual data events. 
Systematic uncertainties due to variations of the track reconstruction efficiency as a function of 
vertex location and also the procedure used to generate the background events are evaluated. 
The individual contributions are added in quadrature to find the final systematic uncertainties 
of 7.3-7.6%. 
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4 Correlation Functions and Near-Side Yields 

The two-dimensional (2D) per-trigger-particle associated yield distribution of charged hadrons 
as a function of | Arj | and | Acp\ is measured for each p^ lg and pj ssoc interval, and in different cen- 

trality classes of PbPb collisions. An example for trigger particles with 3 < p^ is < 3.5GeV/c 
and associated particles with 1 < pip soc < 1.5GeV/c is shown in Fig. [TJ for centralities ranging 
from the 0-5% most central collisions, to the most peripheral (70-80%) events. The 2D correla- 
tions are rich in structure, and evolve with centrality. The and p| ssoc ranges shown in this 
figure were chosen as an example because they demonstrate a good balance of the following 
features. For the most central PbPb collisions, a clear and significant ridge-like structure mostly 
flat in At], and extending to the limit of \Afj\ = 4, is observed at Acp « 0. At mid-peripheral 
events, a pronounced cos(2A</>) component emerges, originating predominantly from elliptic 
flow [10J. Lastly, in the most peripheral collisions, the near-side ridge structure has largely 
diminished, while the away-side back-to-back jet correlations can be clearly seen at Acp n, 
but spread out in At]. A tiny bump at \ At]\ = along the entire \A(p\ range of the correlation 
function is visible for the two most central bins. This is an effect due to low-px "looping tracks" 
mainly at mid-rapidity that pass through the silicon tracker multiple times causing duplicate 
tracks at the same rapidity, but different azimuthal angles. This effect, however, is negligible 
compared to the systematic uncertainties described below. 




Figure 1: Two-dimensional (2D) per-trigger-particle associated yield of charged hadrons as a 
function of \At]\ and \Acp\ for 3 < pf s < 3.5GeV/c and 1 < p| ssoc < 1.5GeV/c, for twelve 
centrality ranges of PbPb collisions at y/s NN = 2.76 TeV. The near-side peak is truncated in the 
two most peripheral distributions to better display the surrounding structure. 
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4 Correlation Functions and Near-Side Yields 



As was done in Ref. 11 J, to quantitatively examine the features of short-range and long-range 
azimuthal correlations, one-dimensional (ID) A<p correlation functions are calculated by aver- 
aging the 2D distributions over a limited region in At] from A to Arj max : 



1 dNP ail 



NMg dA(p Arj r _ 



Mr, 



Al7mi, 



d 2 N? air 



Ntrig dAr] dA(p 



dAt]. 



(5) 



The results of extracting the ID A(p correlations in the short-range (0 < \Ar]\ < 1) and long- 
range (2 < I At] I < 4) regions are shown in Fig. [2] The associatedyield distribution per trigger 

particle is extracted for the same and p| ssoc ranges as in Fig. l] 
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Figure 2: Short-range (0 < | Arj \ < 1, open circles) and long-range (2 < \Arj\ < 4, red closed 
circles) per-trigger-particle associated yields of charged hadrons as a function of \A(p\ for 3 < 
p T 1§ < 3.5GeV/c and 1 < pj ssoc < 1.5GeV/c, for twelve centrality ranges of PbPb collisions 
at y/s NN = 2.76 TeV. The statistical error bars are smaller than the marker size. The systematic 
uncertainties of 7.6% for all data points in the short-range region and 7.3% for all data points 
in the long-range region are not shown in the plots. 



In order to study the short-range A<p correlations in the absence of the flat background in At], the 
ID A(p distribution in the long-range region is subtracted from that in the short-range region. 
The resulting difference of the distributions is shown in Fig. [3] The near-side peak (A(p m 0) 
represents mainly the correlations from jet fragmentation, whereas the away-side region (Acp « 
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7t) is mostly flat and close to zero due to the weak Atj dependence of the away-side jet peak. 
A comparison to the pp data at \/s = 2.76 TeV is also presented, showing a similar structure 
to that in the very peripheral 70-80% PbPb data. However, the magnitude of the near-side 
peak is significantly enhanced in the most central PbPb collisions as compared to pp. Most 
of the systematic uncertainties manifest themselves as an overall change in the scale of the 
correlation functions, with little dependence on Acp and Atj. Therefore, they largely cancel 
when the difference between the short-range and long-range regions is taken. 



CMS L int = 3.9 |ib" 1 PbPb \[s NN = 2.76 TeV 3.0 < p" 8 < 3.5 GeV/c 1 .0 < p* ssoc < 1 .5 GeV/c 




A<|> A<|> A<|) 

Figure 3: The difference between short-range (0 < \Atj\ < 1) and long-range (2 < \At]\ < 4) 
per-trigger-particle associated yields of charged hadrons as a function of \A(p\ for 3 < p™ s < 
3.5 GeV/c and 1 < pj ssoc < 1.5 GeV/c, for twelve centrality ranges of PbPb collisions at y/s NN = 
2.76 TeV. The statistical error bars are smaller than the marker size. The grey bands denote the 
systematic uncertainties. The pp result is superimposed in all panels, for reference purposes. 

The strengths of the near-side peak and away-side region in the \A<p\ distributions from Fig. [3] 
can be quantified by integrating over the two | A^| ranges separated by the minimum position 
of the distribution. This position is chosen as \A(p\ = 1.18, the average of the minima between 
the near side and away side over all centralities. This choice of integration range introduces 
an additional systematic uncertainty in the integrated associated yields. The effect of choosing 
different minima for integration ranges changes the overall yield by an absolute amount of at 
most 0.007. Similar shifts are calculated for each data point, then added as an absolute value in 
quadrature to the other systematic uncertainties. 

Figure|4]shows the integrated associated yields of the near-side peak and away-side regions as a 
function of N par t in PbPb collisions at ^s NN = 2.76 TeV, requiring 3 < Pj 1& < 3.5 GeV/c, for four 
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5 Fourier Decomposition Analysis of the PbPb Data 



different intervals of p| ssoc (1-1.5, 1.5-2, 2-2.5, and 2.5-3 GeV/c). The grey bands represent the 
systematic uncertainties. For easier visual comparison between the most central PbPb results 
and the values one would expect from a trivial extrapolation of the pp results, the latter are also 
represented using horizontal lines covering the full N par t range. The yield of the near-side peak 
increases by almost a factor of two in going from the very peripheral 70-80% to the most central 
0-5% PbPb events, for the lowest pj ssoc interval of 1-1.5 GeV/c. As pip soc increases, the centrality 
dependence of the near-side yield becomes less prominent. This is of particular interest because 
at RHIC energies for p| ssoc down to 2 GeV/c and similar p^ 18 ranges and methodology, although 
for a lower density system (AuAu at y/s NN = 0.2 TeV), there is almost no centrality dependence 
observed [36 J. On both near and away sides, the yield in PbPb matches that in pp for the most 
peripheral events. On the away side, the yield in PbPb decreases with centrality, becoming 
negative for the most central events. This deviation from pp is likely due to a reduction in the 
back-to-back jet correlations for PbPb. 
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Figure 4: The integrated associated yields of the near-side peak (\A(p\ < 1.18) and away-side 
region (\A(p\ > 1.18), requiring 3 < Pj lg < 3.5GeV/c for four different intervals of p2? soc , as 
a function of Np art in PbPb collisions at y/s NN = 2.76 TeV. The statistical error bars are smaller 
than the marker size. The grey bands denote the systematic uncertainties. The lines represent 
the pp results (N par t = 2) superimposed over the full range of N part values. 



5 Fourier Decomposition Analysis of the PbPb Data 

The first Fourier decomposition analysis of long-range dihadron azimuthal correlations for 
PbPb collisions at y/s NN = 2.76 TeV was presented in Ref. d. This analysis was motivated by 
the goal of determining whether the long-range ridge effect was caused by higher-order hy- 
drodynamic flow harmonics induced by the initial geometric fluctuations. This decomposition 
involves fitting the ID A</>-projected distribution for 2 < | At] | < 4 (to avoid the jet peak) with a 
Fourier series given by 

1 HWP air N ( Nmax ) 

where V n & are the Fourier coefficients and N assoc represents the total number of hadron pairs 
per trigger particle for the given \At]\ range and (pj lg , Pj SSOC ) bin. The first five Fourier terms 
(Nmax = 5) are included in both the current fits and those in Ref. [lj. In this paper, the analysis 



5.1 Factorisation of Fourier Coefficients 
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of the Fourier decomposition is extended to the full centrality range, and is performed as a 
function of both pj lg and p| ssoc . 

The Fourier decomposition results have several systematic uncertainties. Because the tracking- 
correction-related systematic uncertainties only change the overall scale of the correlation func- 
tions, instead of the shape, they have only a ±0.8% uncertainty on the extracted Fourier coef- 
ficients (V„a)/ largely independent of n and collision centrality. In addition, the results are 
insensitive to looser or tighter track selections to within ±0.5%. By comparing the Fourier co- 
efficients derived for two different z v tx ranges, |z vtx | < 15 cm and |z vtx | < 5 cm, the systematic 
uncertainties due to the dependence on the vertex position are estimated to be less than ±0.5%. 
Variations from the finite bin width of the Acp histograms contribute the largest systematic un- 
certainty to the analysis, especially for the higher-order components, which are more sensitive 
to the fine structure of the distributions. Reducing the binning of the A(p histograms by factors 
of 2, 4, and 8, the extracted Fourier coefficients vary by ±0.3-2.2%. The effect of including ad- 
ditional higher-order Fourier terms in the fit using Eq. [6] results in changes of at most ±1.0% 
(for n = 5), with Fourier terms up to n = 10 included (N max = 10). The values of additional 
higher-order Fourier terms included in the fit are all consistent with zero. 

Table [2] summarises the different sources of uncertainty for the first five Fourier coefficients. 
These uncertainties are added in quadrature to obtain the total systematic uncertainties, also 
given in Table [2] 

Table 2: Systematic uncertainties of the Fourier coefficients (V n /s) for the first five terms. 



Source Systematic uncertainty of V n & (%) 

Tracking efficiency 0.8 

Vertex dependence 0.5 

Track selection dependence 0.5 

Finite bin width in Acp 0.3-2.2 

Number of terms included in the fit 0.0-1.0 

Total 1.1-2.6 



The fitted Fourier coefficients (V m a) up to n = 5, for two representative low-p T 8 ranges of 
1 < pj ig < 1.5GeV/cand3 < pj ig < 3.5 GeV/c, with p| ssoc fixed at 1-1.5 GeV/c, are presented in 
Fig. [5]for various centrality ranges. The values of V n & peak at n = 2 and then drop dramatically 
toward larger n values at all centralities, although this behaviour is less pronounced for the 
0-5% centrality bin. The error bars show the statistical uncertainties only, while the systematic 
uncertainties are indicated in Table |2] 

5.1 Factorisation of Fourier Coefficients 

If the observed azimuthal dihadron correlations at large A77 are driven only by the single- 
particle azimuthal anisotropy with respect to a particular direction in the event, the extracted 
Fourier coefficients (V„a) from long-range azimuthal dihadron correlations can be factorised 
into a product of the single-particle azimuthal anisotropy harmonics, v„, via 

KA(PT ig 'PT S0C ) = MPT g ) >< Mpf^), (7) 

where v n (pj 18 ) and v n (pj ssoc ) are the harmonics for the trigger and associated particles f20"|, 
respectively. One source of v„ is the collective-flow harmonics arising from hydrodynamic 
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Figure 5: Fourier coefficients Via through V$&, extracted from the long-range (2 < |A^| < 4) 
azimuthal dihadron correlations, for 1 < pj lg < 1.5GeV/c (closed circles) and 3 < pj lg < 
3.5GeV/c (open circles) with p| ssoc fixed at 1-1.5 GeV/c, for twelve centrality intervals. Most 
of the statistical error bars are smaller than the marker size. The systematic uncertainties (not 
shown in the plots) are indicated in Table |2| 
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expansion of the medium (e.g., anisotropic elliptic flow contribution to v%) |37| , particularly in 
the low-px regime where hadron production in heavy ion collisions is thought to be mainly 
from the bulk medium |38| . In addition, for very high pj particles that are predominantly 
produced by the fragmentation of energetic jets, v n could also be induced by the path-length 
dependence of the jet-quenching effect inside the medium ||39l - l44| . This path difference can 
lead to a stronger suppression of the high-px hadron yield along the long axis of the elliptically- 
shaped system than along its short axis, resulting in an azimuthal anisotropy characterised by 
the V2 harmonic. Both scenarios satisfy the factorisation relation of Eq. [7] 

This relation (Eq. [7j| is a necessary ingredient for the extraction of single-particle azimuthal 
anisot-ropy harmonics using the dihadron correlation data, since a Fourier series can be used 
to decompose any functional form by construction. The relation can be tested by taking the 
ratio of V n &(pj lg , p T SS0C ) directly extracted as a function of p^ 18 and p T SS0C (left-hand side of 

Eq. [7)1 to the product of v n (pj 18 ) and v n (pf soc ) (right-hand side of Eq. [tJi, where v n (pj lg ) is 
derived using the fixed, lowest-p T ssoc range (so that the non-flow effect from back-to-back jet 
correlations is minimised) and writing the factorisation relation as 



^trig^ _ VnA{Pj S ,Pj W ) 
Vn{p l T 



The range of pj w is chosen to be 1-1.5 GeV/c. The value of v n (p l j w ) is then calculated as the 
square root of V n ^(p l f w , Pj w ). The ratio of the measured left and right sides of Eq. [7] should 
then be approximately unity if factorisation is valid. 

Figures [6j [7j and [8] show the ratios for n = 2, 3, and 4, over five p T ssoc ranges as a function 

of Pj lg , for both short-range (0 < \At]\ < 1) and long-range (2 < \At]\ < 4) regions. Three 
different centrality intervals 0-5%, 15-20% and 35-40% are presented. The first point of each 
panel equals 1.0 by construction. The error bars correspond to the statistical uncertainties. 

No evidence of factorisation is found for the short-range region (0 < \Atj\ < 1) in any of the 
centrality ranges of the PbPb data. This was expected because of the dominant jet-like corre- 
lations in this region. In contrast, the long-range region (2 < \Atj\ < 4) for mid-peripheral 
15-20% and 35^0% events does show evidence of factorization for V2A to with p| ssoc up to 

3-3.5 GeV/c and pj 18 up to approximately 8 GeV/c. The data are also consistent with factorisa- 
tion for even higher p^ 18 (> 8GeV/c) combined with low pip soc , but the current event sample 

is not large enough to draw a firm conclusion. Note that V n & varies by almost 60% in the p^ 18 
range from 1 to 3.5GeV/c as shown in Fig. |5j whereas factorisation is found to hold to better 
than 5%. This suggests a strong connection between the extracted Fourier coefficients from 
long-range dihadron correlations and the single-particle azimuthal anisotropy harmonics. For 
the most central 0-5% collisions, the ratio for n = 2 deviates significantly from unity, while 
V3A and V4A still show a similar level of factorisation to that of the 15-20% and 35-40% mid- 
peripheral data. This may indicate the existence of other sources of long-range correlations for 
the most central collisions that violate the factorisation relation. Studies of the pp and very 
peripheral (e.g., 60-80%) PbPb data are also carried out. Factorization is found to be invalid, 
as expected, due to the lack of thermalization in these small-size systems [|38| . Higher-order 
Fourier coefficients and a wider pj range can be investigated once larger samples of PbPb data 
are collected. 
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Figure 6: The ratios of Vih{pj , Pj SS0C ) to the product of ^(Px ) anc ^ ??2(pf SS0C ) f° r n = 2 in the 
short-range (0 < \At]\ < 1, open circles) and long-range (2 < \At]\ < 4, closed circles) regions, 
where ^(Pt) is evaluated in a fixed p^, ssoc bin of 1-1.5 GeV/c, for five intervals of p| ssoc and 
centralities of 0-5%, 15-20% and 35^0%. The error bars correspond to statistical uncertainties 
only. 
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centralities of 0-5%, 15-20% and 35^0%. The error bars correspond to statistical uncertainties 
only. 



14 5 Fourier Decomposition Analysis of the PbPb Data 



x 



X 



X 



X 



X 



n=4 CMSL . = 3.9j^b" 1 



1.5 
1.0 

1.5 
1.0 

1.5 
1.0 

1.5 
1.0 

1.5 
1.0 



• 2 < |Ar| 
o < |Ar| 



< 4 o 

< 1 



°oooO° o « 



5°/= 



PbPb \/s^ = 2.76 TeV 



1 .5 < p° aau " < 2^C 



2,0 GeV/c 



° ° 

o°ooo 
2.0 < p* ssoc < 2 T 5 



GeV/c 



o 



oo 



Oq c 



.4- 



2.5 < p* ssoc < 3.0 GeM/c 



■0 < P T < 3.5 



GeWc 



O u I 



4r. ■aysuu j- 
.0 < P T < 5 



GeVc 



^ 



-i , 1 r 1 1 1 r 

15-20% 

°ooo o 



s ' ' I 



o°°°° . 



OO 



o 4 , 



35-40% 



o <> 



o Oo o 







o 



If 



o 



'+1 



•*♦+ 



2 4 6 8 

p*' 9 (GeV/c) 



2 4 6 8 

p tng (GeV/c) 



2 4 6 8 

p tng (GeV/c) 



Figure 8: The ratios of V^(pj, pj SS0C ) to the product of v±(p T ) and v^{p^ soc ) for n = 4 in the 
short-range (0 < |Aj/| < 1, open circles) and long-range (2 < \Ar]\ < 4, closed circles) regions, 
where v^(pj) is evaluated in a fixed p^, ssoc bin of 1-1.5 GeV/c, for five intervals of p| ssoc and 
centralities of 0-5%, 15-20% and 35^0%. The error bars correspond to statistical uncertainties 
only. 
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5.2 Elliptic and Higher-Order Single-Particle Azimuthal Anisotropy Harmonics 



As discussed in Section |5.1[ except for Z7 2 in the very central PbPb events, the factorisation 
relation given by Eq.|7|for long-range (2 < \At]\ < 4) azimuthal dihadron correlations is found 
to be valid for sufficiently low p| ssoc , combined with low p^ 18 , and possibly high py lg as well. 

Therefore, the single-particle azimuthal anisotropy harmonics v n (pj^) can be extracted using 
Eq.|8|withl < p| ssoc < 3GeV/c. Values are found for centralities ranging from 0-5% to 70-80%, 
and presented in Fig. [9] The l-3GeV/c p| ssoc range is chosen in order to reduce the statistical 
uncertainty by utilising as many associated particles as possible over the p| ssoc range where 
factorisation is valid. Data for the most central and most peripheral events are included for 
completeness, although the results for V2 in those events are clearly demonstrated by Fig. [6] to 
be more complicated in nature. The value of V2 is extracted up to pj lg ~ 20GeV/c for all but 

the 2 most peripheral centralities, whereas higher-order v n are truncated at pj lg ~ 10 GeV/ c or 
less for the peripheral data due to statistical limitations. As mentioned previously, factorisation 
is not demonstrated conclusively at very high Pj lg . For the most central 0-5% events, all the 

harmonics are of similar magnitude across the entire pj lg range. The px dependence of all v n 
shows the same trend of a fast rise to a maximum around px ~ 3 GeV/c, followed by a slower 
fall, independent of centrality up to 50-60%. The magnitude of V2 increases when moving 
away from the most central events. At very high Pj lg , sizable signals are observed, which 
exhibit an almost flat pj dependence from 10 to 20 GeV/c for most of the centrality ranges. This 
is not the case for the higher-order harmonics. In order to explicitly investigate the centrality 



dependence of the harmonics, the extracted through v% are also shown in Fig. 10 as a function 



of Npart/ for representative low (1-1.5 GeV/c), intermediate (3-3.5 GeV/c), and high (8-20 GeV/c) 

trig 

p T " ranges. 



A strong centrality dependence of V2 is observed in Figs.|9jand 10 for all p^ 18 ranges, while the 
higher-order harmonics V3-V5 do not vary significantly over the range of N pa rt- This behaviour 
is expected in the context of both the hydrodynamic flow phenomena for lower-px particles |38| 
and the path-length dependence of the parton energy-loss scenario for high-px particles II39I . 
The V2 harmonics are sensitive to the eccentricity of the almond-shaped initial collision re- 
gion that becomes larger for the peripheral events, whereas the higher-order harmonics are 
driven by fluctuations of the initial geometry that have little dependence on the collision cen- 
trality [45 J. In the most peripheral events, the pj dependence of V2 is found to be very different 
from that in the central events, as shown in Fig. |9j In the high-px interval 8-20GeV/c (third 



panel of Fig. 10 1, V2 increases rapidly at low N part (very peripheral). A possible explanation 
is the presence of non-flow effects due to back-to-back jets. Based on the factorised Fourier 
coefficients from long-range dihadron correlations, the single-particle azimuthal anisotropy 
harmonics extracted over a wide range of px and centrality allow a detailed comparison to the- 
oretical calculations of the hydrodynamics and path-length dependence of in-medium parton 
energy loss. 



6 Summary 

The previous CMS analysis of angular correlations between charged particles has been ex- 
panded to cover a wide centrality range of PbPb collisions at y/s NN = 2.76 TeV. As was seen 
previously for central PbPb collisions, the associated yields differ significantly from those ob- 
served in pp interactions. Correlations with both small (0 < \Arj\ < 1) and large (2 < \Ar/\ < 4) 
relative pseudorapidities were again studied as a function of the transverse momentum of the 
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Figure 9: The single-particle azimuthal anisotropy harmonics V2-V5 extracted from the long- 
range (2 < | Arj | < 4) azimuthal dihadron correlations as a function of Pj 1S , combined with 
1 < p| ssoc < 3 GeV/c, for twelve centrality intervals in PbPb collisions at ^/s NN = 2.76 TeV. Most 
of the statistical error bars are smaller than the marker size. The systematic uncertainties (not 
shown in the plots) are indicated in Table |2| 
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Figure 10: The single-particle azimuthal anisotropy harmonics, V2-V5, extracted from the long- 
range (2 < j A?/ 1 < 4) azimuthal dihadron correlations as a function of N par t in PbPb collisions 
at ^/s^ = 2.76 TeV for 1 < pf soc < 3GeV/c in three p^ 8 ranges of 1-1.5, 3-3.5 and 8-20GeV/c. 
Most of the statistical error bars are smaller than the marker size. The systematic uncertainties 
(not shown in the plots) are indicated in Table |2| 



trigger and associated particle pj. The integrated yield of the near-side region shows an in- 
creasing enhancement towards more central PbPb collisions, especially for low-px associated 
particles. 

To further characterise the dependence of the correlations on relative azimuthal angle, a Fourier 
decomposition of the distributions projected onto Acp was performed, as a function of both cen- 
trality and particle pj. Evidence of a factorisation relation was observed between the Fourier 
coefficients ( V n &) from dihadron correlations and single-particle azimuthal anisotropy harmon- 
ics (v„). This holds for pf soc < 3-3.5 GeV/c over the pj lg range up to at least 8 GeV/c in central 
and mid-peripheral PbPb collisions, except for vi in the most central events. The observed 
factorisation is absent in pp and very peripheral PbPb data, indicating a strong connection 
of the observed long-range azimuthal dihadron correlations to the single-particle azimuthal 
anisotropy in heavy ion collisions, such as the one driven by the hydrodynamic expansion of 
the system. The single-particle azimuthal anisotropy harmonics V2 through C5 were extracted 
over a wide range in both pj and collision centrality profiting from the broad solid-angle cover- 
age of the CMS detector. The comprehensive correlation data presented in this paper are very 
useful for studies of the path-length dependence of in-medium parton energy-loss, and pro- 
vide valuable inputs to a variety of theoretical models, including hydrodynamic calculations 
of higher-order Fourier components. 
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